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We consider theoretically the novel technique in magnetic resonance force mi- 
croscopy which is called "oscillating cantilever-driven adiabatic reversals". We 
present analytical and numerical analysis for the stationary cantilever vibrations 
in this technique. For reasonable values of parameters we estimate the resonant fre- 
quency shift as 6Hz per the Bohr magneton. We analyze also the regime of small 
oscillations of the paramagnetic moment near the transversal plane and the frequency 
shift of the damped cantilever vibrations. 
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I. INTRODUCTION 

Magnetic resonance force microscopy (MRFM) based on a cyclic adiabatic inversion (CAI) 
is considered as one of the most promising roads to the ultimate goal of a single-spin detection 
in solids (see, for example, Jl], |J). Typically CAI is generated using the frequency modulation 
of the external radio- frequency (rf) field. In this case, a paramagnetic moment of a sample 
follows the effective magnetic field in the rotating system of coordinates (RSC), and influence 
the cantilever vibrations. 

Recently a new technique called "oscillating cantilever-driven adiabatic reversals" (OS- 
CAR) has been suggested and implemented in ||. In this technique, the cantilever driven 
by an external force causes the CAI of the paramagnetic moment of a sample. The back 
reaction of the paramagnetic moment causes the frequency shift of the cantilever vibrations, 
which is supposed to be detected. The main purpose of this paper is a theoretical analysis 
of the stationary vibrations of the cantilever in the OSCAR technique. Our consideration is 
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FIG. 1: A schematic setup of the system under consideration. Bq + AB is the uniform perma- 
nent magnetic field, B\ is the rotating rf magnetic field, F(t) is an external force acting on the 
cantilever in the z-direction, rhp is the magnetic moment of the ferromagnetic particle, fl is the 
magnetic moment of the paramagnetic cluster, d is the equilibrium distance between the center of 
the ferromagnetic particle and the cluster. 

based on the classical equations of motion for the spin-cantilever system. 

The paper is organized as follows. In Sec. II we introduce the model. The linear OSCAR 
regime is considered in Sec. Ill, and nonlinear regime is analyzed in Sec. IV. A perturbative 
approach and numerical results are presented in Sec. V. In Sec. VI, we analyze the damped 
oscillations of the cantilever in the absence of the external force. In Sec. VII we give a brief 
summary of our results. 



II. HAMILTONIAN AND EQUATIONS OF MOTION 

A schematic setup of the studied system is shown in Fig. 1. A spherical ferromagnetic 
particle with magnetic moment, rap, is attached to the cantilever tip. A small paramagnetic 
cluster with magnetic moment, fl, which must be detected, is placed on the surface of non- 
magnetic sample beneath the tip of the cantilever. The whole system is placed into the high 
permanent magnetic field, Bq + AB, oriented in the positive z-direction. The external force, 
F(t), drives the cantilever vibrations along the z-axis. The transversal rotating magnetic 
field, Bi(t), is applied to the paramagnetic cluster. We place the origin of our coordinate 
system at the equilibrium position of the cantilever tip. 

We consider the cantilever tip as an oscillator with the effective mass, m*, and the effective 
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spring constant, k s . The classical Hamiltonian for the cantilever with the ferromagnetic 
particle and the paramagnetic cluster has the form, 

* =6 + k 4 - ' f ® - s^"- - + aS + « 

where p 2 and 2 are the momentum and coordinate of the cantilever tip, /j,q is the permeability 
of the free space. Putting F(t) = -F cos(z/£) and taking into consideration the finite quality 
factor, Q, of the cantilever we write the equation of motion for the cantilever, 

Z + U c z + ( d + z y +-q Z = fo COS(l/t), (2) 

where o; c = [k s /m*) 1 / 2 is the unperturbed cantilever frequency, / = F /m*, and, 

= 3/i m F 

^ ~ 2vrm* ' [ ' 

Next, we assume that the rf field B ± rotates in the (x,y) plane with the frequency, 

uo = j[B o + B d (0)]. (4) 

Here 7 is the gyromagnetic ratio of the paramagnetic cluster, Bd(z) is the dipole magnetic 
field produced by the ferromagnetic particle at the point of location of the paramagnetic 
cluster, 

and Bd(0) is the value of B d (z) at the equilibrium position of the cantilever, z — 0. The 
equation of motion for the paramagnetic moment jl in the RSC has the form, 

/i = 7 [/Ix 4//]- (6) 

— * 

Here _B e // is the effective magnetic field in the RSC with the x-component B 1 and the z- 
component AB + B' d (z), where B' d (z) is the oscillatory part of the dipole field produced by 
the ferromagnetic particle on the cluster: 

B' d (z) = B d (z)-B d (0)- (7) 
III. THE LINEAR OSCAR REGIME: SMALL OSCILLATIONS OF jl 



In this section we consider the linear OSCAR regime. Suppose that initially an auxiliary 
7r/2-pulse changes the direction of the paramagnetic moment, /3, from +z to +x of the RSC. 
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We also assume that the oscillatory part of the dipole field, B' d (z), is small compared to 
the rf field, B\. Certainly we assume that the unperturbed cantilever frequency, uj c <C 7-Bi, 
to keep the conditions of CAI. In quasi-static approximation a paramagnetic moment, ft, 
follows the effective field B e ff. Putting in (|5|) jl = 0, we obtain for \z\ <C d: 

3fiQm F z( y t) 



fi x {t) « n, fi y (t) = 0, fi g (t) = 



AB 



2nd 4 



(8) 



These equations describe small (linear) oscillations of jl near the x-axis. Substituting the 
last expression in Eq. @ into Eq. (0) we derive an approximate equation for the cantilever 
oscillations, 

z + ufz + -^z = / O cos(i4). (9) 

Here, 

, a a 3/i m F /x / 3{i m F \ f Auj c \ 
u c =u c + Au c , Au c = — — g A5 + , Q = Q 1 + . (10) 

Equation (|9]) describes the motion of the linear oscillator with the effective frequency, 
lu*, and the effective quality factor, Q*. Due to the back reaction of the paramagnetic 
moment on the cantilever the effective frequency and the quality factor of the cantilever 
depend on the permanent magnetic field, AB, (in our approximation, AB <C Bi). If 
AB > —3fj,omF/(87rd 3 ), then both the frequency and the quality factor of the cantilever 
decrease. In the opposite case they increase. 



IV. NONLINEAR ADIABATIC REGIME: ADIABATIC REVERSALS OF p 

To increase the back reaction of jl it is important to provide large oscillations (adiabatic 
reversals) of the paramagnetic moment. In this section we consider stationary vibrations of 
the cantilever in the nonlinear OSCAR regime. It is convenient to write the equations of 
motion in the dimensionless form: 

2 " U + ^ + ^l C * + ^ + ^ 

MJ = (5 - X Z)M y , (11) 
M y ' = eM z -(5- X Z)M x , 
MJ = -eM v , 

where we introduced the dimensionless time r = u> c t; prime means a differentiation over r, 
Z = z/A is the dimensionless coordinate, A = f Q/u^ is the unperturbed (in the absence of 
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the magnetic moment M) amplitude of the stationary cantilever vibrations in the resonant 
regime (when u = u c ), M = p/p is the dimensionless magnetic moment, 5 = 'yAB/uJc. 
The parameter a = A/d is small, a ~ 0.01. The dynamics is controlled by the following 
dimensionless parameters: 

_ 3p m F p 



2vrd 4 QF c 



o 



X= 2™W . (12) 



c 



7#i , . 

£= , p = VUJ c -l. 

Suppose that the paramagnetic moment, M, points initially in the direction of the effec- 
tive magnetic field, B e ff, and the cantilever points in the opposite direction, Z(0) = —1. In 
this case the quasi-static motion of M is given by the expressions, 

M x (t) = 



£ 2 + (5- X Zf 



M y = 0, (13) 
M z (t) = - ■ 

^ £ 2 + {6 _ xZ) 2 

Substituting (0) into the first equation in (|ITD we obtain the nonlinear equation for Z: 

Z" + Z- XXZ +^Z> = ± cos[(l + p)r + U (14) 
U 2 + (xZ) 2 Q Q 



where we neglected the term aZ in the denominator in the third term in the left-hand side 
and put 5 = 0. The third term in Eq. ( |I4D corresponds to the modification of the potential 
energy of the cantilever, due to the interaction with the magnetic moment, M, by the value, 

A 



6u(z) = --^Tixzy. 

Now, we present an approximate "semi-quantitative" analysis of the stationary oscil- 
lations described by Eq. (0). The solution for the stationary driven oscillations of the 
cantilever, described by Eq. (|14D , can be written in the form, 

Z = a(p)sin[(l + p)T + o ]. (15) 

We define the frequency shift which corresponds to the shift of the maximum, a max (pi), of 
the amplitude, a = a(p), caused by the ferromagnetic sample. In order to estimate a max we 
replace, 

Z 2 ~ sin 2 [(l + p)r + # ] = 1{1 - cos[2(l + p)r + 2^ ]} -> 1 
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in the denominator in the third term in Eq. (|i4|), and neglect the term cos [2(1 + p)r + 2$ ] 
because it is non-resonant. Then Eq. (H) takes the form, 



Z" + 



1 



Ax 



e* + \ 



(16) 



The position, p 1; of the maximum of the amplitude, a max , of the driven oscillations (the 
frequency shift) is, 

Ax A 

Pi ~ po = 1 i « — 7=, (17) 

where po — — 1/4Q 2 is the position of the maximum of the amplitude in the absence of the 
paramagnetic sample, and we suppose that x~^> e. 

For estimation of the value of the frequency shift the following parameters where used: 
D — 1.5 x 10~ 7 m is the diameter of the ferromagnetic particle with the volume V = 1.8 X 
1(T 21 m 3 , p m F /V « LIT, k c « 10" 3 N/m, uj c /2tt « 10 5 Hz, A « lnm, d ta lOOnm, 
B\ fa 10 _3 T. For these values of parameters we obtain, 



e w 280, x ~ 2.5 x 10 3 , A w 8.5 x 10^ 5 (p/p B ), a = 0.01, 



(18) 



where p/ps is the paramagnetic moment expressed in units of the Bohr magneton. The 
corresponding frequency shift is, 



Pi ~ Po 



6 x 10- 5 (p/p B ). 



(19) 



This gives the frequency shift -6 Hz per one Bohr magneton. 



V. PERTURBATION APPROACH 



The qualitative estimation presented above can be supported by application of the ap- 
proach based on the perturbation theory developed by Bogoliubov and Mitropolskii in [|J. 
We look for the solution of Eq. (|14|) in the form, 

Z = a(r) cos[y>] + Xui(a, if>), (20) 

where if) = (1 + p)r + i?(t). The function Ui(a, ip) is the sum of the Fourier terms with 
the phases 3ip, 5ip, 7ip, . . .. The amplitudes of these terms decrease with increasing the 
Fourier number, n, as l/(2n + l) 2 . The first non-vanishing term is small and equals to 
u\(a, ip) ~ 0.02 cos (3^). This allows us to neglect the contribution of Ui(a, ip) into the 
expression for Z in Eq. (f20|). 
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The slow varying amplitude, a(r), and the phase, t?(t), in the first order of the pertur- 
bation theory satisfy the two coupled differential equations, 

da A f 2lT x a cos V' sm ipdip a 1 



dr 



27tJ o Je 2 + {xa costp) 2 2 Q Q( 2 + p) 



sin 



d<d 
dr 



1 



8Q 2 



P 



A /" 27r %a cos 2 ^cf^ 



27ra 



1 



(21) 



(22) 



'° yje 2 + (xacos^) 2 aQ(2 + p) 
Note, that the integral in the right-hand side of Eq. ([H]) is equal to zero. The integral in the 
right-hand side of Eq. ( p^) can be expressed through the elliptic integrals as (see Appendix), 
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71 'I 2 x a cos2 ipdip 



4 



k 



E(k)-p 2 kK(k) 



(23) 



le 2 + (xacos-?/>) 2 

where fc = 1/Vl 4- p 2 , K(k) and are the complete elliptic integrals, respectively, of the 
first and second kind, p = e/(ax)- When p 2 <C 1 one can decompose K(k) and E{k) as, 



i2 



E{k) « 1+ C 



1\ A; ,z 



JT(*)«C + (C-1)^ + ... 

where k' 2 = 1 - k 2 ^ p, C = ln(4/fc') m ln(4/p). From Eqs. (f23|) and 
of the integral in Eq. for p <C 1, 



+ . . . 



(24) 



A 



2~ 



Xa cos 2 ^c?^ 
e 2 + (x a cosipy 



2ixa Jo 

Substituting Eq. (p5|) to Eq. (|22"D we obtain 

da a 
aV 



2A 
7ra 



4 



P 



we find the value 



(25) 



dd 
dr 



1 



sg 2 



p 



2A 



7ra 



2Q Q(2 + p) 
4 V p 



sin 



(26) 



aQ(2 + p) 



cos-#. 



We now calculate the position of the maximum of the amplitude, a(p), in the stationary 
regime of driven oscillations using Eq. fl26|), and compare it with the results obtained in 
Sec. |TJ. In the regime of driven oscillations a = const, d = const, and we must solve the 
system of two equations (|26|), where da/dr = and d$/dr = 0. Canceling the phase, we 
have, 



1 



+ g 2 



8Q 



1 2A 

+ P + 



na 



(27) 



a 2 (2 + p) 2 4 

where we neglected the term proportional to p 2 1. The amplitude, a, can be written as 
a = 1 + (3, where (5 <C 1, so that, 



1 



a(2 + p) (l + P)(2 + p) 2 



1-/3 



P 



(28) 
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FIG. 2: The dependence of the amplitude of the driven oscillations of the cantilever on the fre- 



quency detuning, p, obtained using numerical solution of exact equations of motion (11). Solid 
line corresponds to the initial conditions ( |3l[) and the values of the parameters A = 8.5 x 10~ 5 , 
X = 2500, e = 280, a = 0.05. Q = 100, 5 = 0. Dotted line corresponds to the same values of the 

-M z \. Dashed 



parameters but for "inverted" initial conditions [in Eq. (31) M x — ► — M x , M z 
line represents the dependence a(p) with no paramagnetic moment (A = 0). 



Taking the square root from the both sides of Eq. ( p7|) and using Eq. (p8|) we obtain, 

2 



P 



2Q< 



1 2A 



(29) 



2 \8Q 2 7T / 

where we put a ~ 1 in the denominator of the term proportional to A (i.e. we neglected the 
term of the order of (3X). The maximum of the function, (3 = (3(p), can be found from the 
condition d(3{p\)/dp = which yields, 

This is approximately the same value as that given by Eq. (0), obtained from the qualitative 

considerations. The second term in Eq. ( 0O|) describes the influence of the paramagnetic 

moment reversals on the resonance frequency of the cantilever. 

To verify our analytical results we solved numerically the exact equations of motion (PD . 

Fig. 2 (solid line) demonstrates the dependence of the stationary amplitude of the cantilever 

vibrations, a, on the frequency detuning, p. (The stationary amplitude is achieved at r > 

Q.) The initial conditions are taken in the form, 

e ^ x 



Z(0) = -1, Z(0) = 0, M x (0) 



M y (0) =0, M z (0) 



(31) 
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FIG. 3: The dynamics of the projections of the paramagnetic moment, M(r), of the sample with 
the initial conditions (|3l|). The gray line is obtained as a result of the numerical integration of 



Eq. (11), and the black line indicates the quasi-static solution (|13|). For M z (t) both curves almost 



coincide. The parameters are the same as those for the solid line in Fig. 2. 

For these initial conditions at r = the paramagnetic moment, M, points in the direction 
of the effective magnetic field, while the cantilever is displaced in — z-direction ($o = 37r/2) 
from its equilibrium position. 

Fig. 3 demonstrates the motion of the paramagnetic moment, M(r). One can see the close 
correspondence between the analytical and numerical solutions. Note that the frequency 
shift caused by the adiabatic reversals changes its sign if the paramagnetic moment points 
initially in the direction opposite to the effective magnetic field (while Z(0) = —1). Dotted 
line in Fig. 2 depicts this case. We also should note, that decreasing the parameter e (the 
x-component of the effective magnetic field) leads to the violation of the CAI conditions. 
Fig. 4 demonstrates this situation for e = 28. 



VI. DAMPED OSCILLATIONS OF THE CANTILEVER 



The influence of the sample on the cantilever can be measured if one turns off the external 
force acting on the cantilever, and measures the frequency of small damped oscillations of 
the cantilever. In absence of the paramagnetic moment M the frequency of the oscillations 
is independent of time and equals to \Jl — 1/(AQ 2 ). 

We look for the solution of the cantilever vibrations in the form: Z = a(r) cos[r + #(t)]. 
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FIG. 4: The same as in Fig. 3 but for e = 28. 

Then the dynamical equations for the slow varying amplitude, a(r), and phase, i?(t), in the 
presence of the sample and in the absence of the external force take the form, 

M 1 

For p <C 1 Eq. ( p3|) can be written as, 

A? J_ 

The last term in the right-hand side of Eq. (R3) describes a change of the frequency of small 
oscillations of the cantilever caused by the adiabatic reversals of M. From Eq. (|32|) we have 
a(r) = a(0) exp[— r/ (2Q)]. One can see from Eq. that for the initial conditions (|31| ) the 
influence of M results in decrease of the frequency of small oscillations of the cantilever in 
comparison with the case A = 0. For small p (the value of p = ej (ax) ~ exp(t/2Q) increases 
with time) the frequency of oscillations decreases when time increases, as shown in Fig. 5, 
while in the absence of the sample this frequency remains independent of time. We should 
note that in the studied approximation the sample does not influence the amplitude of the 
cantilever oscillations. 



da a 

A f 2n x a cos2 ipdip 



2ira Jo 



\A 2 + ix a cos ipy* 



(32) 
(33) 



2A 
ira 



4 I p . 



(34) 
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FIG. 5: The frequency of small damped oscillations of the cantilever for the initial conditions 



(31) as a function of time, r, for three values of the quality factor, Q. Solid lines are obtained 
using Eqs. (32) and (|34|). The results of exact numerical solution are plotted by the filled circles, 
A = 8.5 x 1(T 5 , x = 2500, e = 280, 5 = 0, a = 0.05. 



VII. SUMMARY 



We have studied theoretically and numerically the stationary cantilever vibrations in the 
novel OSCAR MRFM technique. Our results are based on the application of the classical 
theory for the motion of the cantilever and the paramagnetic moment of a cluster on the 
surface of the sample. We have estimated the resonant frequency shift for the cantilever 
vibrations. For the reasonable values of parameters our estimate is about 6 Hz per Bohr 
magneton. The sign of the shift depends on the initial direction of the paramagnetic moment 
relative to the initial position of the cantilever. We supported our estimation by the analyti- 
cal analysis based on the perturbation theory and by the numerical solution of the equations 
of motion. Our perturbative approach is based on the fact that the influence of the para- 
magnetic moment on the sample is weak (A <C 1). We considered also the regime of small 
oscillations of the paramagnetic moment near the transversal plane (linear OSCAR regime). 
Finally, we analyzed the damped oscillations of the cantilever (without the external force). 
We have shown that the frequency of the damped oscillations becomes time-dependent due 
to the adiabatic reversals of the paramagnetic moment. 
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Appendix 

Here we express the integral in Eq. ( f22|) in terms of complete elliptic integrals, 

r 71 / 2 x acos2 ipdip r n / 2 cos 2 ipdi/j r 71 / 2 (1 — sin 2 ip)dip 

Jo ^fe 2 + {xa cos^) 2 h VP 2 + cos 2 $ Jo ^2 + 1 _ sin 2 ^ 

«r/2 (p 2 + 1) (l - S in 2 ^ - p 2 
^0 VFTT^l-^sin 2 ^ 

where we introduced the notation p — ej {ax)- Splitting this integral in two parts we obtain 
the right-hand side of Eq. (f23|) . 
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